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Integrals necessary for the determination of transition moment matrix elements from experimental 
data have been evaluated numerically by use of vibrational wave functions derived from KKK potentials. 
A power series expansion for the electronic transition moment has been" assumed. The significant 
quantities which can be related to an arbitrary center of expansion are vibrational overlap integrals 
and quantities of the form Ji//,r"i/i,-"r/r. Experimental band intensities and relative populations for 
vibrational levels of the initial electronic state are needed to determine the expansion coefficients. 
Transition moment integrals have been calculated for first ionization from the ground electronic state 
of CO and for the A 2 \\,— X' z ^ + transition of CO + . Comparison of these integrals with previous calcula- 
tions based on Morse functions has shown them to be rather sensitive to the wave-functions [potentials] 
used. Characteristics generally attributed to the r-centroid and related integrals are examined, and some 
limitations on the use of the r-centroid approximation are discussed, following a review of assumptions 
made in the use of that approximation. 
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1, Introduction 

Quantitative determination of transition probabilities 
for diatomic molecules is generally based on the 
"dipole approximation", the bases and limitations of 
which have been stated by Bethe and Salpeter [1]. 
The length form for the transition moment matrix 
element is commonly used for simplicity. Herzberg 
[2| has given the basis for the approximate separation 
of nuclear and electronic multipole moments so that 
only the electronic contribution need be considered. 
The Born-Oppenheimer approximation is assumed 
valid, i.e., 



where r is the internuclear distance and x refers to 
the electronic coordinates. Vibration-rotation inter- 
action and isotope effects are not considered here. 

Intensities of electronic-vibrational bands or radia- 
tive lifetimes of vibrational energy levels are most 
often used to provide the measures of transition 
probabilities. This discussion emphasizes the use of 
intensities. 

The 'integrated intensity of an electronic-vibrational 
band (electric dipole transition) is given by [2] 



*K>taiU, r) ^\ff e (x, r)i// r (r) 



(1) 



Ic>v" = N v >A v > v »hccr v > v »= const. N r '0- A v , v ,,\R v > v » 
for emission, and by 

lv'v"= Const. N v »(T v 'v»\Rv>v»\ 2 



(2) 



(3) 



*A preliminary version of these results was present 
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for absorption. N v > and N& are the populations of levels 
rented at the Moiecukr Spectroscopy v ' and v" ', A r > v » is the transition probability for spon- 
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taneous emission, a v > v » is the wave number of the 
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transition, and rv ri " is the multipole moment matrix 
element whose square is the band strength, S r 'v". The 
transition probability is proportional to the square 
of the matrix element [2] which is given by 

«,v= JMr)*e(r)Mr)<*r (4) 

where the electronic transition moment is given by 



Re(r) = J>r'(*, r)M(x)^e"(x, r)dx. 



(5) 



\fj V ' and i/v are the vibrational wave functions for 
levels v and v'\ i/v and i// P " are the corresponding 
electronic wave functions, and M(x) is the electric 
multipole moment, a function only of electronic 
coordinates (e.g., electronic dipole moment is eXxi). 

For an electric quadrupole or magnetic multipole 
transition, the more complicated expressions analo- 
gous to (5) still reduce to some function of internuclear 
distance, Rt(r) . 

There exists ambiguity about some quantities which 
appear in (2) and (3). A characteristic wave number 
for the band is not clearly defined. If absolute intensi- 
ties are considered then the constants include statis- 
tical weights of the upper and lower electronic states. 
Transition moment and band strength have been 
defined in the literature in several different ways. A 
critical discussion of these and related topics has 
recently been given by Tatum [3]. 

R e (r) usually cannot be calculated from first prin- 
ciples because of insufficient knowledge about elec- 
tronic wave functions, so we assume that R ( > can be 
approximated by 



/? e (r) =#o + #i (r-r () )+/? 2 (r- r ) 2 + . 



(6) 



where r is some meaningful center of expansion for 
the electronic transition. A possible choice for r is 
(rj, + r")/2 which will generally lie within a region where 
both t//, ' and i/v are nonvanishing; however, this choice 
is arbitrary. An r () "characteristic" of the 0-0 transi- 
tion has been suggested in the literature, but this is 
not always a useful choice when \r' € — r n e \ is large, since 
the 0—0 transition may not be easily observable. 
It is tempting to consider (6) a Taylor series and 
to relate the expansion coefficients, /?/, to derivatives 
of R e (r); however, the use of generally imprecise data 
on intensities and vibrational populations, together 
with the inclusion of only a few terms in the expansion, 
makes this procedure uncertain. The transition 
moment can be simply approximated as 



R e (r)=Ro + R { r + R 2 r 2 + . . . 



(6a) 



and the empirical coefficients Rj obtained by using 
the observed values of /, cr, and N in (2) or (3), and 
employing the method of least squares. 

Substitution of (6a) into (4) gives the matrix element 

Rw=Ro(v\i/')+Ri(v\rtr)+R2(v',rV') + ! . . (7) 



the square of which is the band strength. In (7) integrals 
of the form fifjvr'tyvdr are indicated by (i/, rV). 
In the r-centroid approximation [4], the square of (7) 
is approximated as 



S*, |/U-| 2 -«»(?, 



(8) 



The band strength is written as the Franck-Condon 
factor [square of the vibration overlap integral] times 
the square of the transition moment as a function of 
r-centroid which is 



_ (v ', rv") < 
"~ {v\v") ' 



(9) 



Separability of the total wave functions in the Born- 
Oppenheimer approximation does not necessarily 
imply the factorability of the band strength into an 
electronic part times a vibrational part. That factoriza- 
tion is based on an assumed form [r-dependence] for 
the intermediate parameter, the transition moment. 
For constant R e (r), S v > v » includes only a single term, 
Rlq V ' V ", and the Franck-Condon factor contains all the 
relative band-to-band variation. For linear 7? e (r), 
(7) becomes 



or 



R*+=Ro(v' 9 tn+Ri(v',rin 



R**=(v',if)(Ro + Rir*+) 



(7a) 



(7b) 



and Sv'v» has the form of (8) [5, 4]. In the r-centroid ap- 
proximation the assumption is made that the band 
strength is also factorable as in (8) for [slowly varying] 
nonlinear R e (r). This is based on the assumed near 
equality of ratios [4] 



= (v\rv") _ (t/, rV') _ (y\ rV') 



(v',v") (ir\ rv") (i/,r"-V) 



(10) 



up to perhaps n — 10. A consequence of this assump- 
tion is [4] that, for any reasonable function, /(r), 



(v', f(r)v")=f (?,r) (v',v") 



(11) 



which, in particular, always implies (8). Our numerical 
calculations on the ratios of (10) for the ionization 
transition (X 2 S + ) CO+-(X'2 + ) CO and the A-X 
transition in CO + suggest that for nonlinear R e , the 
band strength may not always be factored according 
to (8) [6]. 

An implicit assumption made in this type of calcula- 
tion is that R e must be a slowly varying function of 
r [2]. How slowly varying this should be is not clear. 
Nevertheless, in several cases [7], significant variation 
has been deduced from use of experimental intensities 
and (2) and (8). It is, of course, possible for an appreci- 
able apparent variation in R e over a small range of r to 
arise from inaccurate intensity measurements, inac- 
curate values of q V 'if, or the possibly invalid approxima- 
tion (8). 
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Tyte [8] has cautioned against the casual use of 
experimental intensities in determining the form of 
the transition moment, for the relative intensity distri- 
bution may depend on experimental conditions. There 
are also built-in complications in the use of intensities 
because of blended rotational structure and partial 
overlap of bands. 

The use of radiative lifetimes for determining R?(r) is 
even more complicated than use of intensities as can be 
seen from the following equation 



T V ' = 



1 



1 



v" 



const. 2 cr 3 |/W| 2 

/; v'v" 



(12) 



2. Transition Moment Integrals for X 2 X + (CO + ) 

-X 1 2 + (CO)andA 2 n / (CO + )-X 2 S + (CO + ) 

Our original objective was the tabulation of Franck- 
Condon factors and r-centroids for ionization transi- 
tions of CO and electronic transitions of CO+. Our 
results for X(CO + )-X(CO) and A(CO + )-X(CO) 
showed numerous instances where, within a given v' 
or v" progression, the r-centroids did not vary mono- 
tonically (tables 1—2). Later we became aware of a 
number of cases where this nonmonotonic behavior 
had occurred [9, 7b] for calculations based on both 
Morse and RKR potentials. Meanwhile, there had ap- 
peared a theoretical deduction by Nicholls [10] 
that fW must vary monotonically with respect to 
wavenumber within a progression. 

Nicholls used the Morse function with the Fraser- 
Jarmain approximation of mean exponential coefficient 
[11]. His conclusion, however, does not depend on 
the potential used, but only on the assumption that 
(10) and (11) are valid. In this case/(r) is the potential 
function, and (11) leads to the graphical method for 
determining r> r "[4, 12] 



G{v')-G{v")=U'(f v , v ,,)-U'{-rv>v») 



(13) 



where G is the vibrational term value and U is the 
potential energy. Since (13) always gives monotonic 
values within a progression, it may be necessary to 
revise r-centroids and transition moments determined 
in this manner. 

The transition moment integrals in (7) were evaluated 
by use of vibrational wave functions based on RKR 
potentials [13]. These wave functions have been ob- 
tained by a method previously discussed by Benesch 
et al. [14a]. 

Since the raw wave function is generated as a com- 
bination of two unnormalized segments propagated 
toward each other from either end of the range of 
interest, an accommodation must be made at the 
meeting point. The possible methods of arranging 
for such an accommodation have been discussed else- 
where at length [14]. In the present work, special care 



Table 1. Transition moment integrals for the 

X*2 + (CO+j-X 1 ^ (CO) transition 

Entries have been arbitrarily terminated beyond v" = 10. The 
last two columns are in A. The values of q V 'jf tabulated here are to 
be compared with those based on Morse functions [M. Halmann 
and I. Laulicht, J. Chem. Phys. 43, 1503-9 (1965)]. The integrals 
in (7) are obtained by multiplying the ratios of integrals listed and 
using (v' , v") = qv'v". 









(v' , rv") 


(v' , r-v") 








(v , rv") 


(v' , f z v") 


v' 


v" 


q*ir 


(v', v") 


(v', rv") 






q,'r 


W,v") 


(#', rv") 








0.9644E-00 


1.126 


1.127 












1 





.3566E-01 


1.304 


1.281 


1 


6 


0.1921E-05 


0.666 


0.572 


2 





.9802E-04 


1.839 


1 .654 


2 


6 


.7703E-04 


0.770 


0.735 


8 





.7982E-06 


1.177 


1 . 1 76 


3 


6 


.1414E-02 


0.812 


0.772 












4 


6 


.1776E-01 


0.884 


0.852 





1 


.3383E-01 


0.950 


0.918 


5 


6 


. 1 339E-00 


0.964 


0.929 


1 


1 


.8967E-00 


1.135 


1.138 


6 


6 


.6424E-W) 


1.184 


1.198 


2 


1 


.6923E-01 


1.312 


1.290 


7 


6 


.2038E oo 


1 .356 


1 .338 


3 


1 


.2435E-03 


1.908 


1.691 


8 


6 


.6476E-03 


2.480 


1.926 


4 


1 


.2566E-05 


1.078 


0.739 


9 


6 


.1555E-03 


1.299 


1.149 


8 


1 


.1476E-05 


1.171 


1 . 1 70 


10 


6 


.8092E-06 


2.343 


1.986 


10 


1 


.6424K-06 


1.184 


1.184 


11 


6 


.1117E-05 


1.288 


1.236 


n 


1 


.1251E-05 


1.186 


1.186 


12 


6 


.1308E-05 


1.138 


1.156 


12 


1 


.5070E-06 


1.197 


1.197 






















2 


7 


.7916E-05 


0.728 


0.689 





2 


.1586E-02 


0.865 


0.829 


3 


7 


.1560E-03 


0.760 


0.717 


1 


2 


.6286E-01 


0.954 


0.922 


4 


7 


.2254E-02 


0.813 


0.772 


2 


2 


.8341E-00' 


1.144 


1.149 


5 


7 


.2319E-01 


0.889 


0.860 


3 


2 


.l(X)6K-00 


1.320 


1.299 


6 


7 


.1412E-00 


0.964 


0.927 


4 


2 


.4286E-03 


1.958 


1.718 


7 


7 


.6098E-00 


1.195 


1.211 


5 


2 


.6665E-05 


1.135 


0.830 


8 


7 


.2179E-00 


1.369 


1.350 


6 


2 


.5311E-06 


1.456 


1.471 


9 


7 


.6308E-03 


2.643 


1.979 












10 


7 


.2262E-03 


1.350 


1.206 





3 


.8571E-04 


0.797 


0.754 


11 


7 


.1448E-05 


2.420 


2.008 


1 


3 


.4451E-02 


0.869 


0.834 


13 


7 


.1127E-05 


1.089 


1.120 


2 


3 


.8769F-01 


0.958 


0.926 












3 


3 


.7772E-00 


1.153 


1.160 


1 


8 


.5427E-06 


1.066 


1.101 


4 


3 


.1301E-00 


1.329 


1.308 


3 


8 


.1620E-04 


0.693 


0.622 


5 


3 


.5757E-03 


2.042 


1.758 


4 


8 


.3081 E-03 


0.781 


0.751 


6 


3 


.2249E-04 


1.181 


0.958 


5 


8 


.3286E-02 


0.817 


0.777 


7 


3 


.6232E-06 


1.600 


1.591 


6 


8 


.2738E-01 


0.891 


0.862 


8 


3 


.2231E-05 


1.166 


1.172 


7 


8 


.1445E-00 


0.962 


0.923 


9 


3 


.5205E-06 


1.159 


1.161 


8 


8 


.5837E-00 


1.205 


1.225 


11 


3 


.5928E-06 


1.216 


1.212 


9 


8 


.2392E-00 


1.379 


1.361 












10 


8 


.4082E-03 


3.138 


2.095 




1 


4 
4 


.5492E-05 
.3297K-03 


0.746 
0.807 


0.699 
0.768 


11 


8 


.5033E-03 


1 .350 


1.242 


2 


4 


.8268E-02 


0.874 


0.840 


3 


9 


.2779E-05 


0.741 


0.734 


3 


4 


.1073E-00 


0.961 


0.928 


4 


9 


.3635E-04 


0.721 


0.675 


4 


4 


.7262E-00 


1.163 


1.173 


5 


9 


.4595E-03 


0.764 


0.720 


5 


4 


.1570E-00 


1.337 


1.318 


6 


9 


.4433E-02 


0.817 


0.775 


6 


4 


.7474E-03 


2.103 


1.788 


7 


9 


.3314E-01 


0.898 


0.872 


7 


4 


.4340E-04 


1.226 


1.025 


8 


9 


.1505E-00 


0.962 


0.920 


10 


4 


.5835E-06 


1.146 


1.152 


9 


9 


.5651E-00 


1.218 


1.239 


11 


4 


.9991E-06 


1.146 


1.147 


10 


9 


.2468E-00 


1.392 


1.374 












11 


9 


.2017E-03 


4.116 


2.246 





5 


.5930E-06 


0.785 


0.785 


12 


9 


.5838E-03 


1.420 


1.306 


1 


5 


.2570E-04 


0.752 


0.707 


13 


9 


.1432E-05 


0.601 


6.320 


2 


5 


.7541E-03 


0.806 


0.764 












3 


5 


.1301E-01 


0.880 


0.848 


4 


10 


.6832E-05 


0.770 


0.780 


4 


5 


.1228E-00 


0.963 


0.930 


5 


10 


.6254E-04 


0.716 


0.666 


5 


5 


.6811E-00 


1.173 


1.185 


6 


10 


.6689E-03 


0.761 


0.714 


6 


5 


.1823E-00 


1.346 


1.328 


7 


10 


.5730E-02 


0.818 


0.774 


7 


5 


.7760E-03 


2.235 


1.842 


8 


10 


.3944E-01 


0.906 


0.884 


8 


5 


.5198E-04 


1.304 


1.086 


9 


10 


.1463E-00 


0.956 


0.908 


9 


5 


.1238E-05 


1.870 


1.768 


10 


10 


.5538E-00 


1.229 


1.253 


10 


5 


.5950E-06 


1.115 


1.153 


11 


10 


.2578 K 00 


1.406 


1.388 


11 


5 


.1120E-05 


1.181 


1.192 


13 


III 


.8231 E-03 


1.466 


1.359 



has been taken to insure that the location of the meet- 
ing point be such as (a) to prevent either segment from 
slipping into instability and (b) to provide a relatively 
large value of the wave function ordinate as a basis 
for the rescaling of the separate segments prior to 
joining. Experience has shown that in order to cover 
all cases of mating, with slopes and ordinates large 
and small, it is helpful to introduce a monitoring rou- 
tine for the immediate detection and rectification of 
pathological cases. 
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Table 2. Transition moment integrals for the A 2 !!— X 2 S + transition of CO + 
The last three columns are in A. The integrals in (7) are obtained by multiplying the ratios of integrals listed and using (i/\ v") 
The only deletion in this Deslandres array is 0—13 (q < 0.5 X 10~ fi ). 





■' 


q#* 


(v', rv") 


(v', rV') 


(i/, rV) 


v 


W, v") 


(»', rv") 


(ir\iV) 








0.4237E-01 


1.178 


1.179 


1.181 


1 





.1131E-00 


1.161 


1.162 


1.163 


2 





.1666E-00 


1.145 


1.146 


1.146 


3 





.1797E-00 


1.130 


1.130 


1.130 


4 





.1585E-00 


1.115 


1.115 


1.115 


5 





.1221E-00 


1.101 


1.101 


1.101 


6 





.8519E-01 


1.088 


1.087 


1.087 


7 





.5496E-01 


1.075 


1.074 


1.073 


8 





.3358E-01 


1.063 


1.062 


1.061 


9 





.1959E-01 


1.051 


1.050 


1.048 


10 





.1114E-01 


1.040 


1.038 


1.036 


11 





.6158E-02 


1.029 


1.027 


1.025 





1 


.1518E-00 


1.202 


1.203 


1.204 


1 


1 


.1926E-00 


1.184 


1.184 


1.184 


2 


1 


.9885E-01 


1.167 


1.166 


1.165 


3 


1 


.1422E-01 


1.147 


1.143 


1.139 


4 


1 


.4447E-02 


1.146 


1.153 


1.161 


5 


1 


.4119E-01 


1.125 


1.127 


1.129 


6 


1 


.7890E-01 


1.110 


1.111 


1.112 


7 


1 


.9512E-01 


1.096 


1.097 


1.097 


8 


1 


.9066E-01 


1.084 


1.083 


1.083 


9 


1 


.7431E-01 


1.071 


1.071 


1.070 


10 


1 


.5541E-01 


1.059 


1.059 


1.058 


11 


1 


.3833E-01 


1.048 


1.047 


1.046 





2 


.2505E-00 


1.227 


1.228 


1.228 


1 


2 


.8117E-01 


1.207 


1.206 


1.204 


2 


2 


.2743E-02 


1.205 


1.216 


1.227 


3 


2 


.7248E-01 


1.177 


1.178 


1.179 


4 


2 


.9599E-01 


1.159 


1.159 


1.159 


5 


2 


.4992E-01 


1.143 


1.141 


1.139 


6 


2 


.6767E-02 


1.123 


1.116 


1.109 


7 


2 


.3288E-02 


1.127 


1.136 


1.146 


8 


2 


.2723E-01 


1.107 


1.109 


1.112 


9 


2 


.5362E-01 


1.093 


1.094 


1.095 


10 


2 


.6855E-01 


1.080 


1.080 


1.081 


11 


2 


.6956E-01 


1.068 


1.068 


1.068 





3 


.2522E-00 


1.253 


1.253 


1.253 


1 


3 


.5236E-03 


1.270 


1.296 


1.318 


2 


3 


.1074E-00 


1.217 


1.218 


1.219 


3 


3 


.7054E-01 


1.197 


1.196 


1.194 


4 


3 


.1025E-02 


1.158 


1.137 


1.113 


5 


3 


.3040E-01 


1.170 


1.173 


1.176 


6 


3 


.7028E-01 


1.153 


1.153 


1.154 


7 


3 


.5856E-01 


1.137 


1.136 


1.135 


8 


3 


.2218E-01 


1.121 


1.118 


1.114 


9 


3 


.1022E-02 


1.093 


1.074 


1.054 


10 


3 


.5758E-02 


1.106 


1.113 


1.120 


11 


3 


.2496E-01 


1.090 


1.093 


1.095 





4 


.1736E-00 


1.280 


1.280 


1.280 


1 


4 


.9207E-01 


1.263 


1.265 


1.266 


2 


4 


.8573E-01 


1.240 


1.239 


1.237 


3 


4 


.3553E-02 


1.239 


1.250 


1.261 


1 


4 


.7630E-01 


1.207 


1.208 


1.209 


5 


4 


.5820E-01 


1.188 


1.187 


1.185 


6 


4 


.4001E-02 


1.162 


1.151 


1.139 


7 


4 


.1388E-01 


1.166 


1.171 


1.176 


8 


4 


.4927E-01 


1.147 


1.148 


1.150 


9 


4 


.5485E-01 


1.132 


1.131 


1.131 


10 


4 


.3186E-01 


1.117 


1.115 


1.113 


11 


4 


.7620E-02 


1.101 


1.095 


1.088 





5 


.8628E-01 


1.309 


1.308 


1.307 


1 


5 


.1917E-00 


1.289 


1.289 


1.290 


2 


5 


.4351E-03 


1.218 


1.179 


1.133 


3 


5 


.9682E-01 


1.250 


1.251 


1.252 


4 


5 


.3713E-01 


1.227 


1.224 


1.221 


5 


5 


.7790E-02 


1.224 


1.232 


1.240 


6 


5 


.6208E-01 


1.199 


1.200 


1.201 


7 


5 


.4859E-01 


1.180 


1.179 


1.177 


8 


5 


.5927E-02 


1.157 


1.148 


1.139 


9 


5 


.6402E-02 


1.163 


1.171 


1.179 


10 


5 


.3430E-01 


1.142 


1.144 


1.146 


11 


5 


.4762E-01 


1.128 


1.128 


1.128 





6 


.3206E-01 


1.339 


1.337 


1.336 


1 


6 


.1766E-O0 


1.317 


1.317 


1.316 


2 


6 


.7181E-01 


1.299 


1.301 


1.304 


3 


6 


.6324E-01 


1.273 


1.271 


1.269 


4 


6 


.1862E-01 


1.264 


1.269 


1.274 


5 


6 


.7878E-01 


1.238 


1.238 


1.238 


6 


6 


.1521E-01 


1.214 


1.209 


1.203 


7 


6 


.1157E-01 


1.212 


1.219 


1.225 


8 


6 


.5276E-01 


1.191 


1.192 


1.193 


9 


6 


.4151E-01 


1.173 


1.171 


1.170 


10 


6 


.7227E-02 


1.151 


1.144 


1.135 


11 


6 


.2664E-02 


1.163 


1.176 


1.188 





1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 


1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 


1 
2 
3 
4 
5 
6 
7 
8 
9 

10 

11 



1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 



1 

2 
3 

4 
5 
6 

7 
8 
9 

10 
11 



9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 
9 

10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 
10 



12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
12 
12 

13 . 

13 

13 

13 

13 

13 

13 

13 

13 

13 

13 

13 



q,-i" 



.9157E-02 
.1001E-00 
.1736E-00 
.8726E-03 
.9872E-01 
.6614E-02 
.4152E-01 
.5698E-01 
.5268E-02 
.1402E-01 
.4586E-01 
.3641E-01 

.2012E-02 
.3885E-01 
.1628E-00 
.9119E-01 
.3093E-01 
.5219E-01 
.5274E-01 
.1508E-02 
.5031E-01 
.3841 E-01 
.1301E-02 
.1551E-01 

.3258E-03 
.1104E-01 
.8872E-01 
.1743E-00 
.1598E-01 
.8310E-01 
.4099E-02 
.6895E-01 
.1547E-01 
.1258E-01 
.4873E-01 
.2553E-01 

.3799E-04 
.2314E-02 
.3206E-01 
.1410E-00 
.1289E-00 
.3217E-02 
.8506E-01 
.1063E-01 
.3627E-01 
.4829E-01 
.1165E-02 
.2329E-01 

.2098E-05 
.3642E-03 
.8085E-02 
.6676E-01 
.1717E-00 
.5958E-01 
.4027E-01 
.4330E-01 
.4729E-01 
.3530E-02 
.5099E-01 
.2368E-01 

.9359E-06 
.2899E-04 
.1464E-02 
.2066E-01 
.1093E-00 
.1651E-00 
.1034E-01 
.7871E-01 
.5876E-02 
.6467E-01 
.6190E-02 
.2388E-01 



(v\ rv") 
W,v") 



1.370 
1.347 
1.325 
1.347 
1.283 
1.250 
1.249 
1.227 
1.200 
1.203 
1.184 
1.166 

1.402 
1.378 
1.355 
1.336 
1.304 
1.294 
1.268 
1.283 
1.238 
1.218 
1.182 
1.196 

1.443 
1.411 
1.387 
1.363 
1.352 
1.316 
1.319 
1.279 
1.252 
1.253 
1.229 
1.209 

1.497 
1.445 
1.421 
1.395 
1.373 
1.317 
1.327 
1.292 
1.291 
1.267 
1.218 
1.241 

1.650 
1.479 
1.457 
1.429 
1.404 
1.385 
1.348 
1.339 
1.310 
1.318 
1.278 
1.255 

1.398 
1.575 
1.494 
1.466 
1.438 
1.413 
1.408 
1.360 
1.364 
1.321 
1.284 
1.291 



.1413E-05 
.1713E-03 
.4279E-02 
.4146E-01 
.1489E-00 
.1264E-00 
.1659E-02 
.8401E-01 
.3539E-02 
.4699E-01 
.3228E-01 



1.721 
1.546 
1.506 
1.476 
1.447 
1.424 
1.332 
1.371 
1.318 
1.332 
1.305 



(v\ rV) 

(t/, rv") 



1.368 
1.346 
1.326 
1.373 
1.283 
1.240 
1.252 
1.226 
1.190 
1.209 
1.185 
1.165 

1.400 
1.377 
1.354 
1.338 
1.300 
1.297 
1.266 
1.304 
1.240 
1.215 
1.159 
1.202 

1.439 
1.409 
1.385 
1.364 
1.359 
1.315 
1.333 
1.280 
1.246 
1.260 
1.229 
1.206 

1.489 
1.442 
1.419 
1.394 
1.375 
1.299 
1.328 
1.284 
1.295 
1.265 
1.191 
1.245 

1.609 
1.477 
1.454 
1.428 
1.404 
1.389 
1.344 
1.343 
1.307 
1.333 
1.279 
1.250 

1.421 
1.560 
1.491 
1.463 
1.437 
1.414 
1.418 
1.359 
1.376 
1.321 
1.272 
1.296 



1.675 
1.541 
1.503 
1.474 
1.447 
1.426 
1.300 
1.372 
1.300 
1.335 
1.301 



(i/, rV') 
(i/, rV) 



1.366 
1.345 
1.327 
1.395 
1.284 
1.229 
1.255 
1.225 
1.179 
1.215 
1.186 
1.163 

1.398 
1.375 
1.354 
1.340 
1.296 
1.300 
1.265 
1.323 
1.241 
1.213 
1.133 
1.207 

1.436 
1.407 
1.384 
1.364 
1.365 
1.314 
1.346 
1.280 
1.240 
1.267 
1.230 
1.203 

1.483 
1.440 
1.417 
1.394 
1.376 
1.280 
1.329 
1.275 
1.298 
1.264 
1.160 
1.249 

1.581 
1.475 
1.452 
1.426 
1.404 
1.392 
1.341 
1.346 
1.305 
1.347 
1.280 
1.246 

1.441 
1.548 
1.488 
1.461 
1.436 
1.415 
1.428 
1.358 
1.388 
1.321 
1.260 
1.300 



1.643 
1.536 
1.500 
1.472 
1.446 
1.428 
1.263 
1.373 
1.279 
1.337 
1.298 
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A number of other checking and diagnostic pro- 
cedures were instituted in order to seek out computing 
artifacts and to determine the level of reliability of the 
results. The effects of changes in integration grid size 
were examined. Tables were generated for the point- 
by-point comparison of the function (d^/dr 2 )/^ with 
(U — E) as a measure of compliance with the Schroe- 
dinger equation, while a macroscopic test consisted 
of computing self-overlap integrals and overlaps for 
wave functions of differing vibrational quantum number 
within the same electronic state. In other words, 
tables of fifry'tyvdr have been computed for all the 
values of v' and v" of interest within a single elec- 
tronic state. The self overlap integrals [squared] were 
found to be good to 5 X 10~ 6 and the overlap integrals 
[squared] of orthogonal functions were no larger than 



0.5 X 10~ 6 . Finally, it was determined that calculated 
B v values (i.e., const, (v, 1/r 2 v)) agreed with the cor- 
responding input data values to within 0.0007 to 0.0028 
cm -1 . Accordingly, we rounded the calculated table 
entries to four digits, deleting those for transitions 
having q v < v » < 0.5 X 10" 6 . 

The significant features of tables 1 and 2 are these: 
(a) For numerous vibrational transitions the ratios of 
integrals in (10) are not equal, (b) Within a number of 
v or v" progressions the r V ' V » do not vary monotoni- 
cally. In other progressions the r V ' V » ar e monotonic, 
but are noticeably not a smooth function of v' (or/). 
(c) The onset of nonsmooth behavior correlates with 
minima in the q,-'r"-, whether or not q is small in an 
absolute sense. (Fig. 1, table 3.) 




1 ■ i ■ ■ ■ ■ i ■ ' ■ 
v' = 2 



5 10 

V" 



15 o 




v'=l 



5 10 




FIGURE 1. rv v „ and log q v ^ for A 2 rij-X 2 2\ CO, v'= 1, 2, 3, 4 progressions. 

Calculated values for i/'=14, 15 are based on extrapolation of the potential function for the X state and are not considered reliable. 
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Fraser [4] has stated conditions for the assumed 
validity of the r-centroid approximation: 

(a) fiA(O e — 10 4 ; /jla is the reduced mass of the mole- 
cule in atomic weight units, co ( is in cm" 1 ; 

(b) 0.01 A< |/y-rv»| < 0.25 A; 
(c)rc max is=- 10 in (10); 

(d) q v > V " is not small (though small is not specified). 

Condition (b) is fulfilled for the CO + transition but 
not for the ionization transition, so that in the latter 
case, the r-centroid approximation must be assumed 
to fail under Fraser's criteria. Qualitatively, the non- 
smooth features are more prominent for the ionization 
transition. 

It has been suggested [15] that the functional form 
of the electronic transition moment is more important 
than the choice of wave function in evaluating (4). 



However, the choice of wave function (or, more prop- 
erly, the potential function from which the wave func- 
tion is generated) can also be significant. RKR-based 
q V 'if and r v > v » are compared in table 4 with previous 
values based on Morse functions [16]. This comparison 
does not confirm a statement by Nicholls [17] that the 
parameters (integrals) are insensitive to the exact 
potential used. Flinn [18] and Spindler [19] have also 
shown large differences between RKR and Morse 
integrals, even for small quantum numbers. From all 
these results we conclude that wave functions based 
on RKR potentials should be used where possible so 
that attention can be centered on the form of the transi- 
tion moment. 

Halevi [20] has assumed, as have others before him, 
that most of the contribution to the integrals in (7) 




11 .616-2 .383-1 .696-l\ .250-1 .762-2 .476-1 .266-2 ,364-1 \155-1 .255-1 .233-1 .237-1 .239-1^.323-1 



Table 3. Franck-Condon factors for the A 2 Ilj — X 2 S + system of CO 



Nonmonotonic r-centroids were found, generally, where the ?'s of neighboring bands in a progression differed by a factor of some 10 2 . 
The transition with the smaller q, then, may have an r-centroid which breaks the monotonic series, while the band itself, falling in a hole 
between Condon loci, will be extremely weak or unobserved. Bands which have been observed are underlined. For a compendium of ex- 
perimental data see P. H. Krupenie, "The Band Spectrum of Carbon Monoxide," National Standard Reference Data Series, National Bureau 
of Standards NSRDS-NBS No. 5 (1966). 
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TABLE 4. r/,-,." and r r 'r" for the \-llj — X' 2 1\ C0 + system 

Each entry for grv includes the power of 10 (e.g., -3=10 :t ) by 

which it is multiplied. Most differences between RKR and Morse 
data taken from ref. 16 are small. Only the values showing large 
differences arc compared. 





q,' 


« 


TV* 


(A) 














RKR 


Morse 


RKR 


Morse 


4-1 


4.447-3 


5.796-3 


1.146 


1.136 


2-2 


2.743-3 


3.228-3 


1.205 


1.191 


6-2 


6.767-3 


4.973-3 


1.123 


1.130 


7-2 


3.288-3 


5.155-3 


1.127 


1.116 


1-3 


5.236-4 


5.677-4 


1.270 


1.234 


4-3 


1 .025-3 


5.787-4 


1.158 


1.182 


5-3 


3.040-2 


3.336-2 


1.170 


1.166 


9-3 


1.022-3 


2.34.5-4 


1.093 


1.110 


10-3 


5.758-3 


8.830-3 


1.106 


1.098 


3-4 


3.553-3 


4.186-3 


1.239 


1 .223 


6-4 


4.001-3 


2.730-3 


1.162 


1.174 


2-5 


4.351-4 


3.50.5-4 


1.218 


1.268 


9-5 


6.402-3 


9.00.5-3 


1.163 


1.153 


10-6 


7.227-3 


4.869-3 


1.151 


L.161 


3-7 


8.726-3 


1 .054-3 


1.347 


1 .303 


8-7 


5.268-3 


3.702-3 


1.200 


1.212 


7-8 


1.508-3 


2.289-3 


1.283 


1.253 


10-8 


1.301-3 


4.986-4 


1.182 


1 .205 


.5-10 


3.217-3 


2.976-3 


1.317 


1 .347 


10-10 


1.165-3 


5.736-4 


1.218 


1.252 


1-12 


2.899-5 


4.088-5 


1.575 


1.539 


7-13 


1 .659-3 


3.046-3 


1.332 


1.393 



comes from a narrow region around r,v. This need not 
be so. j\\) v >\\) v »rdr is evaluated numerically as 

V xjj, -'i// r "rAr. Plots of such cumulative sums as a func- 

tion of r (fig. 2) may show regions having large slope 
which indicate large contribution to the integral in 
that r interval. Three general effects occur. Figure 
2(a) describes an integral where there is no contribu- 
tion except very near the average value of r. If the wave 
functions were delta functions then such sums would 
be step functions, and only at a single value of r would 
there be a contribution to the integral. In such in- 
stances the ratios of integrals in (10) would be equal. 
For functions whose product i/vi/v is sharply peaked, 
whether for small or large overlap integral, a plot 
closely similar to figure 2(a) would be expected. 
Figure 2(b) shows smaller oscillations with the most 
significant contribution to the integral coming at large 
r. Figure 2(c)-(g) shows several regions with significant 
slope; the sums show large oscillations about the final 
value of (v' , rv"). 

Figure 2(c)-(g) describes what is likely true of many 
integrals, even for q as large as 0.1 (no integral for 
larger q has been checked, except 0—0, X— X). In 




1.038 1.105 1.172 1.239 1 306 1 373 1087 I 157 1.226 1 246 1.366 1.125 1.192 1.259 1.326 1.393 

r(A) r{ ° M r(A) 

FIGURE 2. ^ i/vi// v »rAr f or A — X, CO + transitions. 

r 
(a) 1-0,(1)) 2-6, <c> 6-2, (d) 7-2, (e) 2-2, (f) 4-3, (g) 2-5. The ordinate is tor smallest r; the ordinate is (i/, rv") for largest r. g,v decreases from (a) = 0.1 to (g) - 0.0004. 



312-248 0-68— 4 
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table 2 and the ratios (10) are not closely the same for 
^2X 10~ 2 . In table 1 the results are erratic. Only 
a detailed examination of all integrals would reveal 
whether the wave functions show delta-function like 
behavior as in figure 2(a). 

Nicholls and Jarmain [12] state: "Comparison 
between the calculated values of r r v' of a band system 
and the ranges of r between the classical turning 
points for each of the levels v' and v" has shown that 
r v > v » lies within, or close to the region which is common 
to both ranges." This is generally true, but James [21] 
has already pointed out examples where f V ' V » lie out- 
side the range of turning points and do not represent 
a meaningful average r. Since products of the wave 
functions do not become insignificantly small till 
0.06 A beyond the turning points, it is not surprising 
if an r-centroid lies beyond the turning points, pro- 
viding tails of both wave functions are involved. Table 
5 shows some of the extreme examples we have ob- 
tained. The A— X, 0-0 and 1-0 transitions illustrate the 
influence of wave function tails. For 0-0 the turning 
points do not overlap by 0.026 A; for 1-0 they overlap 

r" + r' ■ 
slightly. In both instances, r v > v » — -^ — . 

Since the ground states of CO and CO + have nearly 
the same r e and nearly the same potential shapes, 
almost all the ionization transition probability gets 
dumped into a single transition for low v. The 0-0 
transition has q = 0.964. Large values of q are found 
for other members of the Af = sequence, where, for 
example, the 7-7 entry is still above 0.6. 

Table 5. f„ v lying outside turning points for some A — X, CO + 

transitions 
Turning points are taken from ref. 13. fr'r" ana * turning points are 
in A. Each entry for qw includes the power of 10 (e.g., —3= 10~ 3 ) 
by which it is multiplied. 





r,- 1 1-" 


A*n, 


X 2 2 + 






r.i. 


'"ni;ix 


'"miii 


'max 


q,',-» 


0-0 
11-0 
0-6 
0-7 
0-8 
0-9 


1.178 
1.029 
1.339 
1.370 
1.402 
1.443 


1.191 
1.046 

1 1.191 


1.304 
1.634 

1.304 


1.071 
1.071 
0.977 
0.968 
0.961 
0.954 


1.165 
1.165 
1.330 
1.350 
1.370 
1.390 


4.237-2 
6.158-3 
3.206-2 
9.157-3 
2.012-3 
3.258-4 



Halevi [20] has expanded the transition moment in 
a Taylor series centered about the r-centroid and has 
obtained a correction term to R?,(r v > r »)qv>v» which de- 
pends on the second derivative of R e (r v , v „). This cor- 
rection term may have little quantitative significance 
when (10) is not valid. 

3. Conclusions 

We have drawn several conclusions regarding the 
determination of electronic-vibrational transition 
probabilities from a study of two electronic transitions: 

(1) The r-centroid approximation for obtaining 
electronic transition moments and band strengths 



may prove to be not continuously valid along the sig- 
nificant portion of a vibrational progression. The band 
strength may not always be factored as /?f (r r 'r")<7r'r" for 
nonlinear R e . 

(2) Interpretation of the r-centroid as an average r 
is not always meaningful. 

(3) The transition moment integrals are sensitive to 
the potential used; they should be calculated from RKR 
potentials where possible. 

(4) The Fraser criteria are not sufficient to insure the 
validity of the r-centroid approximation. 

(5) The series expansion for the transition moment 
may be used when the r-centroid approximation fails. 



We thank Roger Main for bringing to our attention 
the work by Tatum prior to its publication. The com- 
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to the Computer Science Center of the University of 
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